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ABSTRACT 

The Chandra Deep Field is the region of the sky with the highest concentration of X-ray data available: 4Ms of 
Chandra and 3Ms of XMM-Newton data, allowing excellent quality spectra to be extracted even for faint sources. We 
take advantage of this in order to compile a sample of heavily obscured Active Galactic Nuclei (AGN) using X-ray 
spectroscopy. We select our sample among the 176 brightest XMM-Newton sources, searching for either a) flat X-ray 
spectra (P < 1.4 at the 90% confidence level) suggestive of a reflection dominated continuum or b) an absorption turn- 
over suggestive of a column density higher than ~ 10^"^ cm^^ . Wc find a sample of nine candidate heavily obscured 
sources satisfying the above criteria. Four of these show statistically significant FcKa linos with large equivalent widths 
(three out of four have EW consistent with 1 kcV) suggesting that these arc the most secure Compton-thick AGN 
candidates. Two of these sources are transmission dominated while the other two are most probably refiection dominated 
Compton-thick AGN. Although this sample of four sources is by no means statistically complete, it represents the best 
example of Compton-thick sources found at moderate-to-high redshift with three sources at z=1.2-1.5 and one source at 
z=3.7. Using Spitzersnd i/ersc/ie/ observations, we estimate with good accuracy the X-ray to mid-IR (12/um) luminosity 
ratio of our sources. These are well below the average AGN relation, independently suggesting that these four sources 
are heavily obscured. 

Key words. X-rays: galaxies; Infrared: galaxies 



1. Introduction 

The nature of the X-ray background (XRB) has been con- 
tentious since its discovery 50 years ago (Giacconi et al., 
1962). The Chandra mission confirmed that the background 
is made up of the summed emission from Active Galactic 
Nuclei (Brandt & Hasinger, 2005). The resolved fraction 
of the XRB is about ~90 % in the 0.5-2 keV and 2-5 keV 
bands (Alexander et al., 2003; Xue et al., 2011). Optical 
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spectroscopic follow-up observations show that the peak of 
the redshift distribution of these sources is at z ~ 0.7 — 1 
(Barger et al., 2003; Silverman et al., 2010). At higher ener- 
gies, the limited sensitivity hampers the resolution of a large 
fraction of the XRB. Narrow energy-band source-stacking 
shows that this fraction reduces to -^60% over 5-8 keV, and 
only -50% above 8 keV (Worsley et al., 2005, 2006; Xue et 
al., 2012; Moretti et al, 2012). 

Still, it is at high energies where the bulk of the XRB 

energy density is produced. The peak of the X-ray back- 
ground at 20-30 keV (e.g. Frontera et al., 2007; Churazov 
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ct al, 2007; Moretti et al., 2009) can be reproduced only 
by invoking a significant number of heavily obscured and 
Compton-thick AGN. Compton-thick AGN are those where 
the absorbing column densities exceed 1.5 x lO^** cm~^ and 
thus the attenuation of X-rays by photoelectric absorp- 
tion is enhanced by the scattering on electrons (for re- 
views on Compton-thick AGN properties and surveys see 
Comastri, 2004; Georgantopoulos, 2012). However, the ex- 
act density of heavily obscured AGN required by X-ray 
background synthesis models remains still under dispute 
(Gilli, Comastri & Hasinger, 2007; Sazonov et al., 2008; 
Treister, Urry & Virani , 2009; Ballantyne et al., 2011; 
Akylas et al., 2012). In particular, the intrinsic fraction of 
Compton-thick AGN may vary from 15 to over 35% (see 
discussion in Akylas et al., 2012). Additional evidence of 
a numerous Compton-thick population comes from the di- 
rectly measured space density of black holes in the local 
Universe (see Soltan, 1982). It is foimd that this space den- 
sity is a factor of 1.5-2 higher than predicted by the X-ray 
luminosity function (Marconi et al., 2004; Merloni & Heinz, 
2008), although the exact number depends on the assumed 
efficiency in the conversion of gravitational energy to radi- 
ation. 

The advent of the Integral and SWIFT missions helped 
to constrain the Compton-thick population in the local 
Universe. These missions explored the X-ray sky at energies 
above 10 kcV probably providing the most unbiased sam- 
ples of Compton-thick AGN over the whole sky. Owing to 
the limited imaging capabilities of these missions (carry- 
ing coded-mask detectors), the flux limit probed is bright 
(fi5-55kcv ^ 10^^^ ergcm^^ s^^ ), allowing only the de- 
tection of AGN at very low redsliifts. These ultra-hard sur- 
veys did not detect large numbers of Compton-thick sources 
(e.g. Ajello et al., 2008; Tueller et al., 2008; Paltani et al., 
2008; Winter et al., 2009; Burlon et al., 2011; Malizia et 
al., 2012; Goulding et al., 2011). The fraction of Compton- 
thick AGN in these surveys does not exceed a few percent 
of the total AGN population. In contrast, optical and mid- 
IR surveys yield Compton-thick AGN fractions of between 
10% and 20% (Akylas & Georgantopoulos, 2009; Brightman 
& Nandra, 2011). However, even the ultra-hard surveys 
may miss a fraction of Compton-thick AGN, i.e. those with 
Nh > 2 X 10^'' cm-2 . As Burlon et al. (2011) point out this 
may explain the relative scarcity of Compton-thick AGN 
found in the ultra-hard Integral and SWIFT surveys. 

At higher redshifts, a number of efforts have been made 
to identify Compton-thick AGN. For example, Gilli et al. 
(2010) provide samples of Compton-thick AGN at moder- 
ate redshifts (z ^ 1) using optical spectroscopy and in par- 
ticular the [NeV] emission line. Searches for high-redshift 
Compton-thick AGN in the mid-IR also attracted much 
interest, mainly because the absorbed radiation heats the 
surrounding material and is re-emitted at IR wavelengths. 
The techniques that have been employed include the detec- 
tion of a high 24/xm emission relative to the optical emis- 
sion (e.g. Fiore et al., 2008; Georgantopoulos et al., 2008; 
Fiore et al., 2009; Treister et al., 2009b; Eckart et al., 2010), 
and the presence of a low X-ray-to-mid-IR luminosity ratio 
(see Alexander et al., 2008; Georgantopoulos et al., 2011b; 
Alexander et al. , 2011) Most recently Georgantopoulos 
et al. (2011) proposed the presence of deep silicate fea- 
tures in mid-IR Spitzer spectra, although Goulding et al. 
(2012) argue that the silicate absorption in these sources 
is related to the host galaxies. Nevertheless, the most un- 



ambiguous method for finding Compton-thick sources re- 
lies on X-ray spectroscopy. In X-ray wavelengths, the most 
systematic attempts include these by Tozzi et al. (2006); 
Georgantopoulos et al. (2007, 2009) and more recently 
Brightman & Ueda (2012) all in the Chandra deep fields. 
A number of Compton-thick sources have been individually 
discussed in Norman et al. (2002), Iwasawa ct al. (2005), 
Comastri et al. (2011), Feruglio et al. (2011), GilH et al. 

(2011) Georgantopoulos et al. (2011b) and Iwasawa et al. 

(2012) . In particular, Comastri et al. (2011), provided the 
most direct X-ray spectroscopic evidence yet for the pres- 
ence of Compton-thick nuclei at high redshift, reliably iden- 
tifying two Compton-thick AGN at z=1.536 and z=3.700. 

Here, we attempt to find unambiguous examples of 

Compton-thick AGN at moderate redshifts, extending the 
work of Comastri et al. (2011). The CDFS is one of the 
regions of the sky with the largest accumulation of multi- 
wavelength data available. In particular, it is the area with 
the most sensitive X-ray observations, namely 3 Ms of 
XMM-Newton data and 4Ms of Chandra data, allowing the 
extraction of good quality X-ray spectra even for faint X- 
ray sources such as heavily obscured AGN. Owing to the 
superior XMM-Newton photon collecting power, we first 
select a sample of candidate heavily obscured AGN using 
the 3Ms XMM-Newton observations. For these sources, we 
present a combined XMM-Newton and Chandra spectral 
analysis in order to increase the photon statistics. Our aim 
is to detect unambiguous signs of Compton-thick obscu- 
ration, such as direct detection of a large column density, 
or a flat spectral index, and finally a large equivalent-width 
(hereafter EW) Fe Ka line. Finally, we examine the mid and 
far-IR (up to rest-frame wavelengths of 160/um) properties 
of our sources using data from the Spitzer and Herschel 
missions. The aim is to check whether the IR observations 
independently support a heavy obscuration scenario. We 
adopt Ho = 75kms"^Mpc"\ = 0.3, and Qa = 0.7 
throughout the paper. 



2. Data 

2.1. X-ray 

2.1.1. XMM-Newton 

The CDFS area was surveyed with XMM-Newton during 
different epochs spread over almost nine years. The data 
presented in this paper were obtained by combining the 
observations awarded to our project and observed between 
July 2008 and March 2010 (Comastri et al., 2011) with the 
archival data acquired in the period July 2001 - January 
2002. The total exposure time, after the removal of back- 
ground flares, is w2.82Ms for the two MOS and «2.45Ms 
for the PN detectors. The total area covered is 30x35 ar- 
cmin. The source catalogue contains 339 X-ray sources de- 
tected in the 2-10 keV band with a signiflcancc larger than 
5(7 and a flux density limit of « 6.6 x 10^^"^ erg cm~^ as- 
suming F = 1.7. A supplementary list of 74 sources detected 
with lesser significance is also provided. An extended and 
detailed description of the full data set, including data anal- 
ysis and reduction and the X-ray catalogue will be pub- 
lished in Ranalli et al. (2013). 
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2.1.2. Chandra 

The CDFS 4 Ms observations consist of 53 pointings ob- 
tained in the years 2000 (IMs), 2007 (IMs), and 2010 
(2 Ms). The analysis of the first IMs data is presented in 
Giacconi ct al. (2002) and Alexander ct al. (2003), while 
the analysis of the 23 observations obtained up to 2007 
is presented in Luo et al. (2008). In the present 4Ms sur- 
vey, 740 sources are detected down to a sensitivity limit of 
^ 0.7 X 10~^®crgcm~^s~^ and 1 x 10~^'^ergcm"^ s~^ in 
the hard (2-8 kcV) and soft (0.5-2 keV) band, respectively 
(Xue et al., 2011). The Galactic column density towards 
the CDFS is 0.9 x 10^° cm-^ (Dickey & Lockman, 1990). 

2.2. Infrared 

2.2.1. Spitzer 

The central regions of the CDFS were observed in the mid- 
IR by the Spitzer mission as part of the Great Observatories 
Origins Deep Survey (GOODS). These observations cover 
areas of about 10 x 16.5 arcmin^ using the IRAC (3.6, 4.5, 
5.8, and 8.0 /xm) bands. These data are combined with 
more recent observations of the wider E-CDFS arc^a in the 
SIMPLE survey (Damen ct al., 2011). The combined data- 
set has a 5(7 magnitude limit of [3.6 /im]AB = 23.86, while 
the 3(7 magnitude limit of the central GOODS region is 
[3.6 /im]AB =26.15. The GOODS area in the centre of the 
CDFS has also been imaged with Spitzer-MIPS in the 24 
fim band with a 5(7 flux density limit of 30 /iJy. A much 
wider area, including the entire E-CDFS was imaged as part 
of the FIDEL legacy program (Magnelli et al., 2009) with 
a 5(7 flux density limit of 70 /zJy; we use a combination of 
the two data-sets for this work. 

2.2.2. Herschel 

The far-Infrared data in this work come from the GOODS- 
Herschel survey of Elbaz et al. (2010). This is the deepest 
survey of Herschel using the PACS instrument in both the 
100 and the 160 fim bands, with an integration time of 
more than 15 hr per position. The area covered is a 13x 11 
arcmin field inside the GOODS area. The source detection 
is performed using a 24y^m prior position and the 3ct flux 
density limits are 0.8 and 2.4 mJy in the 100 and 160 yum 
bands respectively. The GOODS-Herschel catalogue we use 
also utilises public data from the HerMES survey (Oliver 
et al., 2012) in the CDFS. We use a 250/Ltm catalogue based 
on 24/um prior positions covering the GOODS-S area, and 
we keep sources with a flux density determination better 
than the 2a limit of 3.5mJy. 



3. Sample definition and X-ray spectral extraction 

3.1. The sample 

We confine our analysis to the brightest sources in the pre- 
liminary XMM-Newton catalogue i.e. those with a detec- 
tion probability of 8ct. There are 194 sources from both the 
main and supplementary lists satisfying the above crite- 
rion. However, for a number of sources we cannot derive an 
XMM-Newton spectrum in any of the cameras, either be- 
cause they are too faint (they lie in areas of elevated back- 
ground or low exposure) or they are confused. The number 



of confused sources is determined using the Chandra imag- 
ing. Therefore, The number of extracted spectra is 176. We 
match the XMM-Newton positions with the Chandra posi- 
tions using a radius of 5 arcsec, and find counterparts for all 
the 176 sources, in the 4 Ms catalogue of Xue et al. (2011), 
or the E-CDFS catalogue of Lehmer et al. (2005). 

We build the multi-wavelength catalogue using all the 
information available and the likelihood ratio method to 
select the counterparts, using the positional uncertainties 
provided in the various catalogues. We first combine the 
SIMPLE catalogue with both the K-selected (Taylor et al., 
2009) and the BVR-sclcctcd (Gawiser et al., 2006) MUSYC 
catalogues (preferring K-sclccted sources in cases where 
they are detected in both catalogues). We find the optical- 
infrared counterparts of the X-ray sources, constraining 
their positions, and then we look for counterparts in the 
FIDEL and 24 /im-prior Herschel catalogues. 

Out of our 176 sources with Chandra counterparts, 136 
have a spectroscopic redshift determination in Szokoly et 
al. (2004), Le Fevre et al. (2005), Norris et al. (2006), 
Ravikumar et al. (2007), Vanzella et al. (2008), Treister, 
Urry & Virani (2009), Balestra et al. (2010), Silverman et 
al. (2010), and Cooper et al. (2011). For 38 of the remain- 
ing sources, photometric redshifts have been compiled from 
Cardamone et al. (2010), Taylor et al. (2009), Rafferty et 
al. (2011), and Luo et al. (2010), Dahlen et al. (2010) while 
two sources have no redshift determination, because they 
are too faint in optical - near-Infrared wavelengths. 

3.2. XMM-Newton spectra 

For each individual XMM-Newton orbit, source counts 
were collected from circular regions with radii between 10 
and 25 arcsec, depending on field crowdedness and source 
brightness, and centered on the 2-10 keV source positions. 
These radii correspond to encircled energy fractions of 
59% and 86%, respectively considering the on-axis XMM- 
Newton PSF at an energy of 4.2 keV (the average energy of 
a source with a spectrum of F = 1). The encircled energy 
fraction does not change abruptly with the off-axis angle 
or the average energy. The above fractions become 55 and 
80% at an off-axis angle of 9 arcmin and for an energy of 
6keV. Local background data were taken from nearby re- 
gions, separately for the PN, MOSl, and M0S2 detectors 
to account for local background variations and chip gaps, 
and avoiding XMM-Newton or Chandra detected sources. 
The areas of the background regions on average have 20 
arcsec radii. The spectral data from individual exposures 
were summed for the source and background, respectively, 
and a background subtraction was made assuming a com- 
mon scaling factor for the source/background geometrical 
areas. Both the PN and the MOS spectra are extracted in 
the 0.5-8 keV energy area. MOSl and M0S2 spectra are 
summed using the FTOOLS MATHPHA task. Response 
and effective area files were computed by averaging the in- 
dividual files using the FTOOLS ADDRMF and ADDARF 
tasks. 

3.3. Chandra spectra 

We used the SPECEXTRACT script in the CIAO v4.2 

software package; to extract the spectra of Chandra sources. 
The extraction radius varies between 2 and 4 arcsec with 
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increasing off-axis angle. At low off-axis angles (<4 arcmin), 
this area encircles 90% of the light at an energy of 1.5 keV. 
The same script extracts response and auxiliary files. The 
addition of the spectral files was performed with the 
FTOOL task MATHPHA. To add the response and auxil- 
iary files, we used the FTOOL ADDRMF, and ADDARF 
tasks respectively, weighting according to the number of 
photons in each spectrum. 

4. X-ray spectral fittings and selection method 

The goal is to identify heavily obscured AGN, via X-ray 
spectral analysis. We use the XSPEC vl2.5 software pack- 
age for the spectral fits (Arnaud, 1996). We employ the 
C-statistic technique (Cash, 1979), which had been specif- 
ically developed to extract spectral information from data 
of low signal-to-noisc ratio. This statistic works on un- 
binncd data, allowing us, in principle, to use the full spec- 
tral resolution of the instruments without degrading it by 
binning. We perform our initial selection in the XMM- 
Newton data both because of its high effective area at high 
energies as well as for its good counting statistics. We fit the 
XMM-Newton spectra using the absorbed power-law model 
PLCABS (Yaqoob, 1997). The advantage of the PLCABS 
model is that it properly takes into account Compton scat- 
tering up to column densities of Nh ~ 5 x lO^'^ cm~^ . 

Our selection method is summarised in the following 
criteria: 

1) The detection of an absorption turnover correspond- 
ing to a column density of Nh > 1-5 x 10^^ cm~^ . For 
a (mildly) Compton-thick source with a column density 
^ 10^* cm~^ , the absorption turnover occurs at rest- 
frame energies somewhat higher than ~8 keV. This implies 
that even with the relatively large effective area of XMM- 
Newton at high energies, it is difficult to detect the absorp- 
tion turnover for marginally Compton-thick sources at low 
redshift. However, at higher redshifts the turnover shifts 
progressively to low energies making the identification of 
Compton-thick AGN more straightforward (Iwasawa et al., 
2012b). For example, for a Compton-thick source with a col- 
umn density of Nh ^ lO^"' cm^^ , the absorption turnover 
would shift to energies about 2 keV at a redshift of z=2, 
an energy region where XMM-Newton has large effective 
area. At a column density of 5 x 10^* cm~^ the turnover 
occurs at an energy of about 20 keV (e.g. Yaqoob, 1997). 
This criterion is reliable only in the case of sources with 
a spectroscopic redshift available. This is because the rest- 
frame energy of the absorption turnover (and hence the 
exact value of the column density) requires knowledge of 
the redshift with high precision. 

2) The detection of a fiat spectral index F < 1.4 at 
a statistically significant level (90% confidence) i.e. the 
90% upper limit of the photon index should not exceed 
F = 1.4. This is an arbitrary selected, albeit extremely 
conservative, cut-off. The average photon index of AGN 
ranges between 1.7 and 2.0 with a standard deviation as 
low as 0.15 (e.g. Nandra & Pounds, 1994; Dadina, 2008; 
Brightman & Nandra, 2011; Ricci et al., 2011; Burlon et 
al., 2011). Therefore, a flat spectrum may be characteristic 
of heavily obscured, reflection-dominated Compton-thick 
AGN (George & Fabian, 1991; Matt et al., 2004). Note how- 
ever that, because of spectral degeneracy, it is likely that 
some sources appear to have flat spectra simply because of 
moderate (~ iQ^^-ss cni"^ ) absorption (see also Corral et 



al., 2011). This degeneracy is pronounced in sources with 
low signal-to-noise spectra. 

3) For the sources which are selected according to the 
above two criteria, we further examine whether the addition 
of a Gaussian component is required by the data. This is 
additional criterion is motivated by the fact that strong 
FcKa lines arc often observed in Compton-thick AGN in 
the local Universe (e.g. Fukazawa et al., 2011), 



5. The Compton-thick candidates 

5.1. XMM-Newton results 

The XMM-Newton spectral fittings yield nine sources with 
either a flat spectrum ^ or a absorption turnover corre- 
sponding to a rest-frame column density greater than Nh ~ 
10^4 cm-2 . The XMM-Newton spectral fits (PLCABS + 
GA in XSPEC notation) of these nine sources are given 
in Table 1. Eight, out of nine candidate Compton-thick 
sources, have a spectroscopic redshift available. One out of 
those (PID-245) has a noisy optical spectrum and a highly 
vmccrtain spectroscopic redshift, (1.864) based on only one 
line (Balestra et al., 2010). For this source there are also 
three photometric redshifts available 2.28, 2.43 and 3.0, 
(Dahlen et al., 2010; Santini et al., 2009; Luo et al, 2010), 
respectively. However, it is possible that the X-ray source 
does not correspond to the counterpart with the available 
optical spectrum (located 1 arcsec away). Iwasawa et al. 
(2012b) detect an FeKalpha line which corresponds to a 
redshift of z=2.68. We chose to use this redshift instead of 
the optical ones. At this X-ray redshift there is no obvious 
optical lin 

For source PID-214, for which there are only photomet- 
ric redshifts available, z=1.5 and z=1.17 (Cardamone et al., 
2010), there is a hint for a line at an energy of E = 2.131q'o3 
keV with EW = 0.3 ±0.1 keV. On the basis of this line, as- 
suming it is related with the FeKa at 6.4 keV rest-frame en- 
ergy, the implied redshift would be z=2.00±0.05. Hereafter, 
we adopt the X-ray redshift for this source. 

At least one source (PID-144) at a (spectroscopic) red- 
shift of z=3.70 is a candidate transmission dominated 
Compton-thick source i.e. is characterised as Compton- 
thick on the basis of an absorption turn-over in its X-ray 
spectrum. This source has first been reported as Compton- 
thick by Norman et al. (2002). A much higher quality X-ray 
spectrum of this source has been reported in Comastri et 
al. (2011). The remaining sources present a flat spectral 
index which may be suggestive of a reflection continuum. 
The FeKa lines provide additional information on the na- 
ture of our sources. Four out of nine sources (PID-66, 144, 
147, 324) present high rest-frame FeKa EW (see table 1). 
These are suggestive of high absorbing column densities 
likely corresponding to Compton-thick AGN. 



^ In addition, there are two more sources which appear to have 

flat spectra but nevertheless the 90% upper limits of their pho- 
ton indices are just above the chosen threshold of F = 1.4. These 
are sources PID-64 and PID-252 at a spectroscopic redshift of 
z=0.516 and 1.893 respectively. 
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Table 1. Initial selection of candidate Compton-thick sources based on XMM-Newton absorbed power-law, (PLCABS), spectral 
fits 



PID 


LID 


XID 


z 


Nh 


r 


EW 


c-stat 


Nh 


Flux 


T obs 


T unobs 

(12) 


Ref 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(13) 


48 


334 


198 


0.298' 




1 9q + U.U4 

J^.^>J_0.09 


0.07+;;;- 


2909/2498 


o q+O.lY 
''••J-0.14 


40.0 


0.68 


0.88 


- 


66 


126 


191 


1.185^ 


1 Q+45.0 


-0.371^:^6 


0.36l»;« 


3015/2498 




1.2 


1.5 


2.2 




144 


265 


412 


3.700^ 


inr + TS.O 


^■■^-1.4 


'-'•^'-'-0.30 


2845/2498 


67.0t^2l:^ 


1.55 


3.7 


40 


a,c,d,e 


147 


176 


257 


1.537^ 


18.6+^-7 


n 4Q+0.59 
^•^'5-0.43 




2881/2498 


61.01^9:° 

/IS+2-5 
48_4.o 


6.3 


2.2 


3.5 


a,b,d 


214 


55 


83 


2.00-^ 


21.7+?-? 


1 Q+0.15 


n 91 +0.09 
'-'•^J^-0.09 


2874/2498 


19 


6.1 


8.0 




222 


428 


686 


0.424^ 


U.OU_Q 29 




<0.27 


2896/2498 


9 n+0.52 
^•*J-0.41 


7.0 


0.3 


0.3 




245 


184 


22 


2.68^ 


45.0l?r.°o 


0.93l°-|? 


<0.62 


2820/2498 


91 n+331 

t»j..u_22.6 


3.0 


1.5 


5.0 


e 


289 


128 


193 


0.607^ 


0.561°:^? 


0.82i°;i| 


11+010 


2850/2498 


q 9+0.4 
'^■''-0.4 


15.9 


1.2 


1.4 




324 


398 


634 


1.222^ 


< 0.63 


0.69l°;f8 


n 71+0.58 
<J- ' J--0.45 


2940/2498 


6-011^ 


3.6 


1.2 


1.2 


f 



The columns arc: (1) XMM-Newton ID. (2) Chandra ID from Luo ct al. (2008). (3) Chandra ID from Xuc ct al. (2011). (41 
redshift; two and three decimal digits denote X-ray and spectroscopic redshift respectively; reference: ' Balestra et al. (2010), 
Szokoly et al. (2004), X-ray redshift from current work. (5) Nn column density in units of 10^^ cm~^ . (6) Photon index. (7) 
FeKa equivalent- width in units of keV; the energy of the line has been fixed at a rest-frame energy of 6.4 keV. (8) c-statistic 
value / degrees of freedom (9) Column density in units of 10^^ cm~^ for a photon-index fixed at P = 1.8. (10) 2-10 keV flux in 
units of 10~i^ ergcm"^ s"' . (ll) 2-10 keV obscured luminosity in units of 10*^ ergs"'^ . (12) Unobscured luminosity (estimated 
by removing the absorption component) in units of 10*^ ergs"'^ . The flux and luminosities are derived having P free. (13) 
previous X-ray reference: (a) Tozzi ct al. (2006); (b) Georgantopoulos et al. (2007); (c) Norman ct al. (2002); (d) Comastri et 
al. (2011); (e) Iwasawa et al. (2012); (f) Georgantopoulos et al. (2011b). All errors refer to the 90% confidence level. 



Table 2. Joint XMM and Chandra spectral fits of the full candidate Compton-thick sample 



Power-law Reflection 



PID 


cts 


cts 


z 


Nh 


r 


EW 


c-stat 


P 


EW 


c-stat 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


48 


7736 


4172 


0.298 


9 i;i+0.23 


1 no+0.U4 
1.00_Q gg 


O.lOiro'a 


3622/3010 


9 c:r+0.03 


<0.013 


4078/3011 


66 


1044 


215 


1.185 




1.8 


U.OU_,) 28 


3614/3010 


1.8 


0.55«:^3^ 


3666/3012 


144 


848 


407 


3.700 


AO 7+44.0 
48.7_27.0 


-'-•^"-0.47 


61 +" ''-^ 


3477/3010 


1 ni+o.i9 
-'-•^-'--0.14 


0.5641;;-;; 


3477/3011 


147 


1207 


593 


1.537 


18.61^7.0 


0.46t-^^ 


0.431;;-;; 


3485/3010 


1 4Q+O.IO 


0.38t-.- 


3487/3011 


214 


2921 


1921 


2.0 


19.8ll« 


o.87l°:li 


0.30l"o:l"o 


3498/3010 


1 09+0.05 


o.3oi°:Sg 


3582.5/3011 


222 


1694 


1021 


0.424 


1 Q+0.40 
^•°-0.40 


1 ri+O.lQ 
^•^J^-0.19 




3477/3010 


1.8 


<0.15 


3702/3012 


245 


460 


259 


2.68 


io6.olJf,:« 


'"-0.82 


<0.37 


3388/3010 


1 cq+0.10 


<0.44 


3404/3011 


289 


3790 


2295 


0.607 


i-3oi;;:^i 


1 90 +0.06 
-'-•^"'-().12 


n is+"'08 

"•-^«-0,()8 


3330/3010 


1.8 


<0.10 


4263/3012 


324 




iOl 


i.222 


< ().()() 


().9()_;;;^;; 


-'-■''-(), r,4 


:512o/:5()l() 


i.8 


o-^i-!;;3H 


:ilo7/:i()i2 



The columns are: (1) XMM-Newton ID. (2) Net XMM-Newton counts (all three modules). (3) Net Chandra counts. (4) redshift; 
two and tree decimal digits denote photometric and spectroscopic redshift respectively. (5) Nh column density in units of 10^^ 
cm~^ . (6) Photon index. (7) Equivalent-width in units of keV; the energy of the line has been fixed at a rest-frame energy 
of 6.4 kcV. (8) c-statistic value / degrees of freedom (9) Photon-index for the PEXRAV model. (10) Equivalent-width of the 
FeKa line in the PEXRAV model. The line energy is fixed to 6.4 keV. (11) c-statistic value/degrees of freedom. All errors 
refer to the 90% confidence level. 



5.2. Joint XMM-Newton /Chandra spectral fits 
5.2.1. Power-law and Fe line model 

To increase the photon statistics, we present the combined 
XMM-Newton and Chandra spectral fits. These are given in 
Table 2. Note that, there are about 70 XMM-Newton and 
Chandra observations spanning a period of more than ten 
years. Then, in principle one should use 70 different nor- 
malisations which is not feasible because of computational 
limitations. Therefore, for the sake of simplicity, the XMM- 
Newton and the Chandra power-law normalizations have 
been tied to the same value. 

In agreement with the XMM-Newton spectral analysis, 
presented in the previous section, four sources (PID-66, 
144, 147 and 324) present FcKq lines with high (rest-frame) 
EW (>0.4 keV). Therefore, the joint Chandra and XMM- 



Newton analysis corroborates that these four sources have 
a good probability for being Compton-thick. In Fig. 1 we 
present the joint Chandra and XMM-Newton FN spectral 
fits. 

5.2.2. Reflection model 

Moreover, we fit a reflection model, Magdziarz & Zdziarski 
(1995) (PEXRAV -I- GA in XSPEC notation). The param- 
eters in the PEXRAV model are: the Fc abundance, the 
inclination angle of the reflecting slab, the cut-off energy of 
the incident power-law and finally the slope of the incident 
power-law. We found that the first three parameters remain 
practically unconstrained owing to the large uncertainties. 
Therefore, we fix the Fe abundance to 1, the inclination 
angle of to 45 degrees while we assume that the incident 
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Fig. 1. The X-ray Chandra (circles) and XMM-Newton PN (squares) X-ray spectra of the candidate Conipton-thick sources PID- 
66, PID-144 and PID-147 and PID-324. The first and second row give the transmission and reflection dominated AGN respectively 
according to the combined Chandra and XMM-Newton spectral fits. Note the high spectral response of XMM-Newton compared 
to Chandra at high energies. 



power-law spectrum has no high energy cut-off (see e.g. 
Dadina, 2008). 

For four sources, the derived photon index is either 
very flat (PID-66, T = 0.80 ± 0.2) or very steep (PID-222, 
r = 2.95 ± 0.08 (PID-289, T = 2.90 ± 0.06), (PID-324, 
r = 2.93^";36). As these values are many a discrepant with 
the intrinsic spectral slopes encountered in AGN (see e.g 
Nandra & Pounds, 1994; Dadina, 2008), we choose to fix 
the photon index to F = 1.8. The results are shown again 
in Table 2. 

5.3. More complex spectral models for the four 'secure' 
Compton-thick AGN candidates 

For the four 'secure' candidate Compton-thick sources, we 
present some additional spectral models for the joint XMM- 
Newton and Chandra spectral fits. First, we are exploring 
the exact energy of the FeKa line. This is motivated by 
the results of Iwasawa et al. (2009) who find ionized Fe 
emission in many Ultra-luminous IRAS galaxies. Moreover, 
we investigate the effect of leaving the ratio between the 
Chandra and the XMM-Newton power-law normalization 
free. Note that we keep the Fe line normalization constant, 
as this is believed to arise far away from the black hole 
in type-2 AGN (Shu et al., 2011). Therefore, the present 



model has two more free parameters: the line energy and 
the Chandra power-law normalization which is not tied up 
to the XMM-Newton power-law normalization (Table 3). 
Note that the normalization of the Chandra and XMM- 
Newton power-laws are consistent within the errors (see 
discussion in section 5.2.1) 

A further model contains an additional power-law whose 

normalization has been set to 3% of the primary power-law. 
The slope of this power-law has been fixed to the same value 
as the primary (transmitted) power-law. This component 
models the scattered (unabsorbed) component along the 
line of sight which is often detected in type-2 AGN (e.g. 
Turner et al., 1997). In this model again the energy of the 
FeKa line is left as a free parameter. The same holds for the 
normalization of the Chandra power-law which is not tied to 
the XMM-Newton normalization. The results are presented 
in Table 4. The spectral fits are in general consistent with 
the simple power-law plus FeKa line fits. The exception is 
source 66 where the photon-index becomes flatter at the 
expense of a lower column density. In this case, because of 
the flatter power-law, the EW becomes higher reaching a 
value of 1 keV in the case of Chandra . In good agreement 
with the power-law model presented in section 5.2.1, the 
four sources (PID-66, 144, 147 and 324) present FeKa line 
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Table 3. Joint XMM-Newton and Chandra spectral fits to the four 'secure' Compton-thick sources Power-law + FeKa line 
{PLCABS+GA) leaving the relative XMM-Newton and Chandra normalization and the energy of the line free. 



PID Nh F E EWx EWc c 

(1) (2) (3) (4) (5) (6) (7) 

"1)6 78.0ls-o 1.52i;;-'^;^ 6.45lo4p 0.23to OSTfJ^^^s 3560/3008 

144 80.0l^-^ l-79l°-;° 6.49+°-^^ 0.47l°'^° 0.56l°-^^ 3475/3008 

147 2Q.Qts°Q 0.57l°J| 6.40l°;lg 0.35l°;?° 0.43^°:^^ 3477/3008 

324 <0.45 0.76lg;l4 6.38lg;^^ O-GQlHa l-15l?:go 3402/3008 



The columns are: (1) XMM-Newton ID. (2) column density (in units of 10^^) (3) photon index as derived from the PLCABS 
model. (4) FeKa line energy (5) XMM-Newton rest-frame EW of the FeKa line. (6) Chandra rest-frame EW of the FeKa line. 
(7) c-statistic value and degrees of freedom. Note: the model is the same as that used in Table 2 with the difference that the 
FeKa line energy is free and also the Chandra power-law normalization is not tied to the XMM-Newton one. 



Table 4. Joint XMM-Newton and Chandra spectral fits to the four 'secure' candidate Compton-thick candidates, using a power-law 
-I- FeKa line -I- scattered emission model (PLCABS-|-GA-|-PO), leaving the relative XMM-Newton and Chandra normalization 
and the central energy of the line free 



PID 


Nh 


r 


E 


EWx 


EWc 


C 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


66 
144 
147 
324 


28.3+?-J 

iqr + lS'.O 

8.011° 
<0.42 


r, r,+0.24 

^■^-0.32 

1-001^3? 
0.72ig;» 


6 45+"-^'^ 

O-*O-0.32 

6.571°:?° 
6.39l°:°« 


^^•^^^-0.20 
U.DO_o.3o 

0.28l°;i 

"•O-'--0.25 


gq+u «o 

U.»»_Q 4Q 

0.70l°;fo 
0.30i°i° 

1 99+0.90 
J-.^^_0.70 


3551/3007 
3458/3007 
3490/3007 
3402/3007 



The columns are: (1) XMM-Newton ID. (2) column density (in units of 10^^) (3) photon index as derived from the PLCABS 
model. (4) FeKa line energy (5) XMM-Newton rest-frame EW of the FeKa line. (6) Chandra rest-frame EW of the FeKa line. 
(7) c-statistic value and degrees of freedom. 



EW higher than ~0.4 keV corroborating their classification 
as good Compton-thick candidates. 



6. Co-added X-ray spectra 

Although there are at least four sources with large FeKa 
line EW, it is possible that some of the remaining sources 
have large EW which fail detection owing to the limited 
photon statistics. In order to answer this question, we de- 
rive the co-added X-ray spectrum. The objective is to de- 
tect faint FeKa emission that cannot readily be detected 
in a single source (e.g. Iwasawa et al., 2012b). The data 
stacking is a straight sum of the rest-frame spectra of in- 
dividual sources. The individual spectra are binned so that 
they have a 2-10 keV rest-frame energy range with 200 eV 
channel width. The data, after background subtraction, are 
corrected for an efficiency curve determined by the response 
matrix and auxilliary response file. Finally, each energy bin 
is corrected for redshift. 

In Fig. 2, we present the co-added rest-frame X-ray 
spectra for two groups of sources. The first group includes 
the four 'secure' Compton-thick AGN candidates (PID-66, 
144, 147 and 324). The second group includes the remain- 
ing sources (PID-48, 214, 222, 245 and 289). A difficulty 
with the second group is that two of the rcdshifts (PID-214 
and 245) are based on the X-ray spectra and thus may be 
more ambiguous. The first group of sources displays a very 
prominent FeKa line at a rest-frame energy of 6.4 keV. The 
second group possibly shows an emission feature but at a 
higher energy of about 7 keV. If confirmed, this could imply 
the presence of an ionized FeKa line. 



7. IR properties 

7.1. Spectral Energy Distributions 

We construct SEDs (Spectral Energy Distribution) for the 
full sample of our flat spectrum sources, to constrain the 
predominant powering mechanism (AGN or star formation) 
in the mid/far IR part of the spectrum. The SEDs are also 
used to obtain an accurate estimate of the 12/im nuclear 
infrared luminosities and thus the X-ray to mid-IR lumi- 
nosity ratios which are often considered as good diagnostics 
for heavily obscured sources. 

We use data from UV wavelengths to the far-IR (where 
available). Four out of nine sources do not have a far-IR 
measurement available from Herschel. The data have been 
modelled using the code originally developed by Fritz et 
al. (2006). The code has been updated by Feltre et al. 
(2012). The SED fitting is based on a multicomponent anal- 
ysis, e.g. Vignali et al. (2009), Hatziminaogiou et al. (2009) 
Pozzi et al. (2012). The observed UV to far-IR SED is 
de-composed in three distinct components: a) stars, hav- 
ing the bulk of the emission in the optical/near-IR b) hot 
dust, mainly heated by UV/optical emission due to gas ac- 
creting on to the supermassive black hole c) colder dust 
heated by star formation. The stellar component has been 
included using a set of simple stellar population (SSP) spec- 
tra of solar metallicity and ages ranging from «1 Myr to «8 
Gyr. A common value of extinction is applied to stars of all 
ages, and a Calzetti et al. (2009) attenuation law has been 
adopted (Ry = 4.05). In order to model the emission at 
observed wavelengths greater than 24/im, the SED fitting 
also includes a component from colder, diffuse dust, likely 
heated by start formation processes, as discussed below. 
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Fig. 2. The co-added (rest-frame) X-ray spectra for two sets of sources. Left panel: the four 'secure' candidate Compton-thick 
AGN (see text). Right panel: the remaining five sources. 
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Fig. 3. Tlie spectral energy distributions of the nine Compton-thick candidate sources. The red (dot), blue (short dash) and green 
(long dash) lines denote the stellar, torus and starburst component respectively, while the black line denotes the sum of the three 
components. 



This is represented by templates of known starburst galax- 
ies (Vignali et al., 2009). The SED fits are shown in Fig. 3. 
The four 'secure' candidate Compton-thick sources appear 
to have star-forming components. However, the strength of 
the star-forming component in PID-66 and PID-147 is am- 
biguous because of the lack of far-IR data in these sources. 
Note that the use of templates of nearby starbursts to 
model the far-IR emission and hence star-formation rates. 



may present some limitations. This is because there may 
be discrepancies between typical IR SEDs at high redshift 
compared to local sources (e.g. Elbaz et al. , 2011; Nordon 
et al., 2012). 

Next, we derive the exact star- format ion rates and spe- 
cific star-formation rates, i.e. the ratio of star-formation 
over the stellar mass. The far-IR luminosity is probably 
the most reliable tracer of star-forming activity (Kennicutt 
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& Evans , 2012). The IR photons are emitted by the dust 
surrounding young stars, which is heated by their ultra- 
violet radiation. The star-formation rate is derived from the 
total IR luminosity (8 — 1000/im) which is ascribed to star- 
formation according to the SED decomposition. We use the 
relation between the IR luminosity and SFR in Murphy et 
al. (2011): 

SFR[MQyr-i] = 3.88 x 10-^^L(8 - 1000;um)[ergs-^] (1) 

Then, the specific star-formation rate, sSFR, is given by 

logsSFR[Gyr-i] = logSFR[MQyr-V logM*[M0] + 8.77(2) 

The stellar mass, M* is estimated from the optical part of 
the SED (for details sec Rovilos et al., 2012). Table 5 sum- 
marises the IR luminosities, star-formation rate and stellar 
masses for the heavily obscured AGN in our sample. The 
typical errors on the 12/im are about 20% apart from PID- 
48, PID-66, PID-222 and PID-245 where the uncertainties 
can be as high as a factor of two. In most cases, the typical 
errors on the total IR luminosity are between 10 and 20%. 
In the cases where no far-IR data are available, the errors 
become large (factors of 2-3). Finally, for the stellar masses 
the errors are of the order of 10%. 

7.2. X-ray to IR luminosity ratios 

The detection of a low X-ray to mid-IR luminosity ratio has 
been widely used as the main instrument for the detection 
of faint Compton-thick AGN which cannot be easily identi- 
fied in X-ray wavelengths (e.g. Goulding et al., 2011). This 
is because the mid-IR luminosity (e.g. 12/im or 6/im) is a 
good proxy of the AGN power, as it should be dominated 
by very hot dust which is heated by the AGN (e.g. Lutz et 
al., 2004; Maiolino et al., 2007). At these wavelengths the 
contribution by the stellar-light and colder dust heated by 
young stars should be small. Gandhi et al. (2009) present 
high angular resolution mid-IR (12/zrn) observations of the 
nuclei of 42 nearby Seyfert galaxies. These observations pro- 
vide the least contaminated core fluxes of AGN. These au- 
thors find a tight correlation between the near-IR fluxes 
and the intrinsic X-ray luminosity. Although the Spitzer 
observations do not have the spatial resolution to resolve 
the core, the SED decomposition allows us to derive with 
reasonable accuracy the nuclear IR luminosity. In Fig. 4, 
we present the obscured X-ray luminosities against the 12 
fxm luminosities. All four candidate Compton-thick sources 
have low Lx/Li2 ratios. Therefore, it appears that the use 
of the (obscured) X-ray to mid-IR (12/im) ratio is a reli- 
able diagnostic for the presence of heavily obscured sources. 
Note that the unobscured X-ray luminosity of at least two 
sources lies below the Gandhi et al. (2009) relation. This 
may suggest that these are Compton-thick despite the ab- 
sence of high-EW FeKa lines. 

7.3. Star formation 

It has been suggested that highly obscured AGN at X- 
ray wavelengths may be associated with large amounts of 
star-formation (see for example Georgakakis et al., 2003; 
Rovilos et al., 2007; Mainieri et al., 2011). Alexander et al. 
(2005) find that many sub-mm galaxies, which have very 
high rates of star-formation, are associated with heavily 



absorbed sources in X-ray wavelengths. A possible explana- 
tion for such a link would be that the nuclear star-forming 
region is associated with the X-ray absorbing screen. From 
Fig. 3, it appears that in the four 'secure' Compton-thick 
sources, the rest-frame mid-IR wavelengths (around 6/im) 
are dominated by the torus emission. Only two of the four 
'secure' Compton-thick sources have long-wavelength data 
(PID-144, PID-324), and the star- forming component is 
dominated in both cases. In Fig. 5 we plot the specific star- 
formation rate as a function of redshift. The four candidate 
Compton-thick sources have a specific star-formation which 
roughly follows the main star- forming sequence as described 
by Elbaz et al. (2010). Moreover, it appears that the star- 
formation properties of the Compton-thick AGN are not 
different from those of X-ray selected AGN in general (see 
e.g. Mullaney et al., 2011; Santini et al., 2012; Rosario et al., 
2012; Rovilos et al., 2012). Therefore, it appears that there 
is no link between the presence of a Compton-thick nucleus 
and enhanced star-formation activity in our sample. 

8. Discussion 

8.1. Limitations in the sample selection and completeness 

We compile a sample of candidate Compton-thick AGN 
searching for either a) flat spectrum sources i.e. those with 
r < 1.4 (at the 90% confidence level) which is suggestive 
of a reflection dominated spectrum (e.g. George & Fabian, 
1991) or b) those showing an absorption turnover suggest- 
ing a rest-frame column density of > 10^"^ cm~^ . The deep 
flux limits of the present observations facilitates the detec- 
tion of Compton-thick sources. This is because Compton- 
thick AGN are diSicult to detect because they are very 
faint in the X-rays for the same intrinsic luminosity, mean- 
ing that many Compton-thick sources at a given luminosity 
will be missed because of the effective survey flux limits. 

We find nine candidate! sources among the 176 XMM- 
Newton spectra examined. Out of these, four show FeKa 
lines with large equivalent-widths and hence these have 
higher probability of being Compton-thick. Our sample 
is by no means complete. For example, at faint fluxes 
there may be sources which have fiat spectra, albeit with 
larger error bars on F and hence fail our selection criterion. 
Incompleteness may also be introduced by the errors in un- 
certain rcdshifts. This is important in the cases where the 
Compton-thick classification is based on the detection of an 
absorption turnover. 

Next, we discuss whether there are apparently flat X-ray 
spectrum sources which are not genuinely reflection domi- 
nated Compton-thick sources. AGN have a steep power-law 
photon index with a slope of F ^ 1.8 and a dispersion of 
a = 0.15 (e.g. Nandra & Pounds, 1994). Sources which 
have a significantly flatter photon-index have a large prob- 
ability of being Compton-thick. However, there are flat- 
spectrum sources which are probably not Compton-thick. 
This is because there is a degeneracy between the photon- 
index and the column density and therefore a moderate 
column density or a complex/multiple absorber can be mis- 
taken as a flat photon index. For example, in table 1, we 
sec that source PID-66 (z=1.185) shows a photon-index of 
F = -0.37^0:6 and Nh = 12ti^ cm'^ (3015/2498) when 
both the photon-index and the column density are left free. 
When the photon-index is fixed to the more realistic value 
of F = 1.8, the implied column density becomes Nh « 
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Fig. 4. The absorbed X-ray (2-10 keV) X-ray luminosity against the 12/im luminosity. The large circles correspond to the best 
candidate Compton-thick sources as suggested by the large equivalent-width FeKa lines. The typical errors are of the order of 
30% and 20% for the IR and X-ray luminosity respectively, including the uncertainties in the model fitting. The hatched diagram 
represents the la envelope of the local (Gandhi et al., 2009) relation. Open circles correspond to X-ray luminosities corrected for 
absorption (see table 1). For two sources (PID-222 and PID-324) the absorbed luminosity equals the unabsorbed luminosity. The 
dotted line corresponds to a factor of 30 lower X-ray luminosity as is typical in many Compton-thick nuclei. 



100^20 ^ 10^^ cvar"^ , while the difference in the Cash statis- 
tic (3018/2500) is not statistically significant. This degener- 
acy can be seen even in sources with large numbers of pho- 
ton counts. Corral et al. (2011) present an X-ray spectral 
analysis of --300 AGN from the XMM-Newton bright sur- 
vey. There are several sources (e.g. XBSJ134656. 7+580315 
at a redshift of z=:0.373), at bright fluxes (f2-io > 10^^"^ 
ergcm^^s^^ ), and hence with excellent photon statistics 
which present a flat power-law as the best-fit. Corral et al. 
(2011) point out that there are no iron lines with large EW 
detected in these sources. If the photon-index is fixed at 
r = 1.9 the resulting column density is of the order ~ 10^"^ 
cm~^ and therefore these sources, although certainly heav- 
ily obscured, are not Compton-thick. 

8.2. The X-ray spectrum of Compton-thick AGN 

Additional ambiguities in the selection may be introduced 
by uncertainties in the X-ray spectrum of Compton-thick 
AGN. Brightman & Ueda (2012) perform a selection of 
Compton-thick sources among the Chandra sources in the 



4Ms CDF-S data using their own spectral models. These 
take into account Compton-scattering, the geometry of the 
circumnuclear matter as well as the scattered nuclear light. 
They find a number of sources with a large scattering com- 
ponent (up to 44%) and hence steeper spectra. The amount 
of scattered light increases with increasing column density. 
These amounts of scattered light are well above those de- 
tected in local Compton-thick AGN (e.g. Comastri et al., 
2010). Evidently, Compton-thick sources with steep spectra 
would have avoided detection in our selection criterion. 

The strength of the FeKa line introduces an additional 
uncertainty. The presence of an FeKa line with a large EW 
is considered to be the 'smoking gun' for the presence of 
a Compton-thick nucleus. Comastri et al. (2010) present 
Suzaku observations of a few Compton-thick sources. They 
find narrow, FeKa lines (EW^l-2 keV) due to neutral or 
mildly ionized gas in all of them. Fukazawa et al. (2011) 
present Suzaku observations of a sample of nearby Seyfert 
galaxies among which many are Compton-thick. They find 
that all Compton-thick sources present ¥eKa lines with an 
EW exceeding 0.5 keV. However, there have been rare ex- 
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Fig. 5. The specific star-formation rate as a function of rcdshift for eight of the candidate Compton-thick sources. The shaded area 
gives the star-formation main-sequence according to Elbaz et al. (2010). Above and below this area lie the starburst and quiescent 
galaxies respectively. Source 222 is not ploted as there is no available rest-frame photometry above a rest-frame wavelength 24/um. 



Table 5. Infrared properties of heavily obscured sources 



PID 


z 


Ll2 




AGN frac. 


logM* 


SFR 


sSFR 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


48 


0.297 


1.66 


1.05 


0.25 


10.51 


1.26 


-0.48 


66 


1.185 


4.21 


1.40 


0.22 


11.21 


1.39 


-1.05 


144 


3.700 


14.60 


33.7 


0.01 


10.87 


2.77 


0.67 


147 


1.530 


7.37 


3.33 


0.10 


10.81 


1.76 


-0.28 


214 


2.00 


9.5 


6.5 


0.95 


10.52 


1.90 


-0.99 


222 


0.424 


0.11 


_t 




10.64 






245 


2.68 


17.7 


6.6 


0.05 


10.34 


2.05 


0.56 


289 


0.605 


1.25 


2.0 


0.22 


10.81 


1.43 


-0.61 


324 


1.222 


75.3 


8.3 


0.36 


10.84 


1.96 


-0.11 



The columns are: (1) XMM-ID ; (2) rcdshift (3) rest-frame 12/im luminosity in units of solar luminosities (xlO^") ; (4) total 
IR luminosity in units of solar luminosities (xlO^^); (5) AGN fraction of the total IR luminosity (6) logarithm of stellar mass 
in units of solar masses (7) logarithm of star-formation rate in units of solar masses/yr; (8) specific star-formation rate in units 
of Gyr"'^. ^: no photometry above 24/xm rest-frame. 



amples of Compton-thick sources with small EW. One of 
these examples is the Broad-absorption-Line QSO Mrk231. 
Its BeppoSAX spectrum (Braito ct al., 2004) shows that 
it is Compton-thick (with a column density of ~ 10^"* 
cm~^ ), but it has an FeKa line with an EW of only O.SkeV. 
Moreover, recent Suzaku observations showed a decrease 
in the covering fraction of the absorber (Piconcelli et al., 
2012). 

8.3. Previous studies of heavily obscured AGN in the CDFS 

Tozzi et al. (2006) first derived a sample of candidate 
Compton-thick sources in the IMs CDF-S observations 



based on X-ray spectral fits. They derived a sample of 20 
candidate Compton-thick sources. Georgantopoulos et al. 
(2007) performed the same exercise compiling a sample of 
18 sources in the same dataset but using slightly different 
spectral models. Only nine of the sources are common in 
these two samples, demonstrating that the uncertainties in 
this technique are considerable. Prom the candidate heavily 
obscured sources in Tozzi et al. (2006) or Georgantopoulos 
et al. (2007), only two sources are found in our present 
heavily obscured XMM-Newton sample. 

Comastri et al. (2011) used the 3Ms XMM-Newton data 

in the CDF-S attempting to confirm that some of the 
Compton-thick candidates in the above samples are indeed 
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secure Conipton-thick AGN. The excellent quality XMM- 
Newton data and in particular its ability to detect the FeKa 
line allowed the above authors to confirm the best two ex- 
amples of Compton-thick candidates in the high rcdshift 
Universe: PID-144 and PID-147 at redshifts of z=1.53 and 
z=3.7 respectively. Note that PID-144 has been first re- 
ported as Compton-thick by Norman et al. (2002) on the 
basis of IMs Chandra spectroscopy. 

For PID-144, Comastri et al. (2011) fit an absorbed 
power- law model finding a photon-index of F = 1.48;to;4o 
with a column density of 6^2 ^ 10^^ cm~^ i.e. favouring a 
transmission-dominated heavily obscured source. Moreover, 
they detect a strong FeKa line with an EW of 84O+420 
eV. Our joint XMM-Newton /Chandra spectral fit (ta- 
ble 2) suggests a transmission-dominated highly obscured 
spectrum with Nh ~ Slg x 10^^ cm^^ and a photon in- 
dex of F = 1.26tg:^7 fully consistent with the results of 
Comastri et al. (2011). The FeKa line EW is somewhat 
lower (GlOlgQg) eV but fully consistent with the above re- 
sults. 

In the case of PID-147, Comastri et al. (2011) find that 
a good fit to the data is provided by an absorbed power- 
law model with a column density of ^ 10^^ cm~^ and a 
photon index of F = —0.11 ± 0.22; the iron line has a very 
large EW (1870^45^5 eV). Our combined XMM-Newton and 
Chandra spectral fits (Table 2) give a higher column den- 
sity of I9I7 ^ -'-0^^ cm~^ but with a steeper photon-index 
r w 0.46j^o!20' from which a smaller equivalent width of 
430^200 results. Further exploration of the spectra sug- 
gests that the difference can be ascribed to degeneracies 
in the parameter space, coupled with the different treat- 
ment of statistic (C-stat vs. x^) and binning of the spectra, 
which may result in different local minima for the statistic 
parameter. 

Brightman & Ueda (2012) report 20 secure Compton- 
thick AGN. Out of these, five have been reported as 
Compton-thick in Tozzi et al. (2006) while three in 
Georgantopoulos et al. (2007). The two Compton-thick 
sources reported by Comastri et al. (2011) (144 and 147) 
and in our paper have column densities just below 10^** 
cm~^ in the spectral fits of Brightman & Ueda (2012). 
Recently, Iwasawa et al. (2012) performed a search for heav- 
ily obscured AGN in the same XMM-Newton sample with 
this presented here. They are using an X-ray colour-colour 
diagram based on the rest-frame 3-5 keV, 5-9 keV and 9- 
20 keV bands. They find a sample of seven high-redshift 
(z> 1.7) candidate heavily obscured sources having rest- 
frame 9-20 keV excess emission. Some of them are good 
Compton-thick candidates either on the basis of a high-EW 
(>1 keV) FeKa line (e.g. PID-114 at z=1.806) or directly 
on the basis of large column densities, Nh ~ 10^^, (PID- 
245, 252). 

8.4. Comparison with X-ray background synthesis models 

Next, we compare with the predictions of X-ray back- 
ground synthesis models (e.g Gilli et al., 2007; Treister, 
Urry & Virani , 2009; Ballantyne et al., 2011; Akylas et al., 
2012) In our sample, we find four candidate Compton-thick 

sources two of which appear to be transmission-dominated. 
Comparing with the models of Gilli et al. (2007)^^, we find 

^ www.bo.astro.it/~gilli/counts.html 



that the number of Compton-thick sources is ^7 in our 
field down to the flux limit of the present survey /(2-10) ~ 
1 X 10~^® ergcm~^s~^ . The model of Treister, Urry & 
Virani (2009)'^ predicts a number of 8 Compton-thick in the 
same area. Finally, the model of Akylas et al. (2012) yields 
about 11 Compton-thick sources in our field. Note that all 
the above estimates refer only to transmission dominated 
(logNn = 24 - 25) Compton-thick AGN. If we include the 
reflection-dominated AGN (logNn = 25 — 26) in the model 
of Gilli et al. (2007) the total number of Compton-thick 
AGN rises to « 11. Interestingly, the model of Akylas et al. 
(2012) has a lower intrinsic fraction of Compton-thick AGN 
compared to that of Gilli et al. (2007): 15% compared to 
33%. However, the higher number of Compton-thick AGN 
predicted by Akylas et al. (2012) at the relatively bright 
fluxes probed here, is a consequence of the stronger re- 
flection component predicted by this model. The fact that 
the number of the Compton-thick candidates found here 
is systematically lower than those predicted by all X-ray 
background synthesis models may suggest that our sample 
suff'ers from incompleteness. 

8.5. Concluding remarks 

The selection of Compton-thick AGN on the basis of ei- 
ther an absorption turnover (transmission dominated) or 
a flat spectrum (reflection dominated) appears to be sufii- 
ciently robust. The presence of a high EW FeKa line can 
be considered as the final diagnostic for the presence of a 
Compton-thick AGN despite some possible exceptions men- 
tioned above (e.g. Mrk231). On the basis of the above diag- 
nostic we isolate four sources as excellent candidates for be- 
ing Compton-thick. These are probably the best Compton- 
thick candidates at moderate to high redshifts found so far. 
Still, the number found here is a factor of at least two lower 
compared to the predictions of X-ray background synthesis 
models. This may imply that a number of Compton-thick 
AGN (particularly those at the faintest fluxes) fail to be 
classified as Compton-thick. It is likely that these lie among 
our remaining five flat spectrum sources. Then we did not 
classify them as 'secure' Compton-thick because we failed 
to detect strong EW FeKa lines. The fact that a few sources 
lie below the Gandhi et al. (2009) (unabsorbed) X-ray to IR 
luminosity relation may point towards this scenario. This 
stresses the need for complementary IR methods together 
with X-ray spectroscopy, in order to better understand the 
properties of Compton-thick AGN. 

9. Summary 

We report on an X-ray spectral study of the 176 brightest 
sources in the XMM-Newton survey in the CDFS. The aim 
is to identify very highly obscured (Compton-thick) AGN. 
Our methodology consists of looking for sources which have 
either a) an absorption cut-off characteristic of a high col- 
umn density or b) a flat spectrum (the 90% upper limit of 
the photon index should be lower than 1.4). After the selec- 
tion of our candidate sources, we are additionally looking 
for the presence of a strong FeKa line which is considered to 
be the trademark of a highly obscured source. Our results 
can be briefly summarised as follows: 

^ agn.astroudoc.cl 
* indra.astro.no.gr 
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— The XMM-Newton spectra select nine candidate heavily 
obscured sources. We separate a group of four sources 
which have large iron line EW and therefore these have 
higher probability of being Compton-thick. Adding the 
Chandra data to the spectral fits corroborates our previ- 
ous results. Two of the sources arc most likely transmis- 
sion dominated (PID-66 and PID-144) while the other 
two are reflection dominated (PID-147 and PID-324). 

— Although our sample is by no means statistically com- 
plete, it represents the most well studied examples of 
Compton-thick AGN at moderate to high redshifts. In 
particular, the redshifts of the four most secure can- 
didate Compton-thick sources are z=1.185, 1.222, 1.53 
and 3.7. One of our candidate Compton-thick sources 
(PID-66) is presented here for the first time. The XMM- 
Newton only spectra of two sources (PID-144 and PID- 
147) have been reported in Comastri et al. (2011) while 
their combined XMM/Chandra spectrum is presented 
here for the first time. Finally, one source (PID-324) 
has been reported in detail in Georgantopoulos et al. 
(2011). 

— The X-ray to mid-IR (12/im) luminosity ratio of the four 
candidate Compton-thick are well below the average re- 
lation of Gandhi et al. (2009), independently suggesting 
heavy obscuration. 

— Spitzer and Herschel observations of our four candi- 
date Compton-thick sources derive star-formation rates 
between about 25 and 1000 Mq yr"^. Their specific 
star-formation rates are consistent with those of nor- 
mal galaxies suggesting that heavy obscuration is not 
related to enhanced star-formation. 
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